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Abstract 

The cosmological moduli problem for relatively heavy moduli fields is rein- 
vestigated. For this purpose we examine the decay of a modulus field at a 
quantitative level. The modulus dominantly decays into gauge bosons and 
gauginos, provided that the couplings among them are not suppressed in the 
gauge kinetic function. Remarkably the modulus decay into a gravitino pair is 
unsuppressed generically, with a typical branching ratio of order 0.01. Such a 
large gravitino yield after the modulus decay causes cosmological difficulties. 
The constraint from the big-bang nucleosynthesis pushes up the gravitino 
mass above 10 5 GeV. Furthermore to avoid the over-abundance of the sta- 
ble neutralino lightest superparticles (LSPs), the gravitino must weigh more 
than about 10 6 GeV for the wino-like LSP, and even more for other neutralino 
LSPs. This poses a stringent constraint on model building of low-energy su- 
persymmetry. 
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Moduli stabilization is one of the long standing problems in the efforts to connect super- 
string theory to the real world. Recent development of the flux compactifications [f] implies 
that some of the moduli as well as the dilaton are stabilized at ultra high energy scale close 
to the Planck scale. However others will remain light compared to the Planck mass. Their 
masses are expected to be not far from the electroweak scale in the context of low-energy 
super symmetry. The remaining moduli will play important roles in phenomenology and 
cosmology. 

In fact, as is termed the cosmological moduli problem [2], such remaining moduli would be 
cosmo logical embarrassment. It is likely that the moduli fields have Planck scale amplitudes 
in the early universe and thus their coherent oscillation would dominate the energy density 
of the universe. Since their interactions are very weak, typically suppressed by the Planck 
mass, they are long-lived. If the masses are around the electroweak scale they decay much 
later than 1 sec, releasing huge entropy. This will completely upset the success of the 
primordial nucleosynthesis, one of the most important triumphs of the big-bang cosmology. 

A resolution of the cosmological moduli problem is to invoke a relatively heavy moduli 
with mass of 10 5 GeV or more [3-6]. 1 Then the moduli fields decay before the nucleosynthe- 
sis commences and will not spoil it. It has been recently recognized that such heavy moduli 
masses can naturally be realized when the moduli have superpotential of exponential type, 
which are typically generated by some non-perturbative effects. For instance the superpo- 
tential of one exponential with a constant (the KKLT-type superpotential [11]) gives the 
modulus mass about 8n 2 heavier than the gravitino mass [12-16], whereas the superpotential 
of two exponentials (race-track superpotential [17]) gives an additional factor of 8ti 2 to the 
modulus mass [18]. 

Thus a natural setting will be that the moduli have masses by several orders of magnitude 
larger than the gravitino. The gravitino can again be heavier than the superparticles in 
the minimal supersymmetric standard model (MSSM), depending on how supersymmetry 
breakdown is mediated to the MSSM sector. 

The purpose of this paper is to reinvestigate cosmological issues of the heavy moduli 
fields. To this end, we will examine the decays of a modulus field, in particular, into a 
gravitino pair at a quantitative level. A special attention is paid to the helicity ±1/2 
components of the gravitino. We will show that the decay amplitude is proportional to the 
F-auxiliary field expectation value of the modulus field. For the heavy modulus field this is 
suppressed by the ratio of the gravitino mass to the modulus mass. However, the net result 
does not suffer from any suppression factor in general. A typical branching ratio of the 
modulus decay into gravitinos is at a % level, which is natural from the counting of degrees 
of freedom. A previous estimate used in Refs. [19-21] is not applicable for a general case. 

The non-suppressed branching ratio of the decay into the gravitinos has striking impacts 
to cosmology. First of all, given a large gravitino yield the gravitino decay would spoil the 



Other ways out have also been proposed. The moduli fields may be fixed at a enhanced symmetry 
point so that the initial amplitude of the moduli oscillation may be small [7]. Late time entropy 
production [8,9] may take place to dilute the unwanted relics. It has been suggested recently that 
the moduli fields may decay rapidly in the thermal bath [10]. Finally one may seek for string 
compactifications with all moduli being stabilized at a scale close to the string scale. 
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success of the big-bang nucleosynthesis. This reincarnation of the gravitino problem [22] 
implies that the gravitino should weigh more than 10 5 GeV. A severer constraint on the 
gravitino mass will be obtained, however, by the over-abundance of the neutralino dark 
matter, provided that one of the neutralinos in the minimal supersymmetric standard model 
(MSSM) is the lightest superparticle of the whole theory and the R-parity is conserved so 
that the LSP is absolutely stable. The requirement that the neutralino LSP abundance 
does not exceed the observation leads the lower bound of the gravitino mass of the order 
10 6 GeV for a weak scale LSP mass when the LSP is wino-like so that the annihilation is 
most effective. On the other hand, the neutralinos produced by the gravitino decay will get 
thermalized if the gravitino mass is larger than 10 7 GeV. 

Such a heavy gravitino is very embarrassing when one attempts to realize weak scale 
superparticle masses in the MSSM sector. We will speculate possible resolutions to this 
problem at the end of the paper. 

The rest of the paper will be devoted to detail the aforementioned results. 

We begin by reviewing some properties of moduli fields. In the following, we assume, for 
simplicity, that there is one modulus field under consideration. As was mentioned previously, 
it is a scalar field in a chiral supermultiplet. We denote it by X. Its properties are governed 
by the Kahler potential K(X, X*), a real function of X and its complex conjugate field X*, 
and the superpotential W{X), a holomorphic function of X. In the following we will use 
the Planck unit where the reduced Planck scale Mp\ ~ 2.4 x 10 18 GeV is set to unity unless 
otherwise stated. The kinetic term of the X field is given by 

Kxx.dpX'&'X (1) 

where the indices X and X* in the Kahler metric Kxx* represent the partial derivatives of 
K with respect to X and X*, respectively. In string models, one often obtains the Kahler 
potential of the type 

K(X,X*) = -n\og(X + X*) (2) 
where n is an (integer) constant. In this case the Kaher metric becomes 



n 

(x + x*y 



K ™ = lv ■ W ( 3 ) 



The F-auxiliary field of the X is given by 

F x = -e K/2 {K X x*)-\W x + K X W)* (4) 

When it is much heavier than the gravitino mass m 3 / 2 = (e K ^ 2 W), the modulus mass is 
given by 

m x = (e K l\Kxx*r l Wxx) = - (^f) , (5) 

where (• • •) stands for a vacuum expectation value (VEV). Supergravity corrections of order 
m 3 / 2 have been neglected. Finally, as is known well, the stationary condition of the potential 
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with respect to the X field implies that the VEV of Fx is generically of order m^ 2 /mx (in 
the Planck unit) when the X field takes a VEV around the Planck scale. Whether it really 
survives non-zero or not will be model dependent. For interesting cases such as the KKLT- 
model and the two race-track model, the VEV of F x field is non-vanishing and indeed of the 
order my 2 /m X - In the following we assume this is the case for the modulus field X under 
consideration. 

We now discuss various decay modes of the modulus field. Let us first consider the decay 
into gauge bosons. The coupling of the modulus field to vector supermultiplets (i.e., gauge 
bosons and gauginos) is conducted by the gauge kinetic function S(X). The relevant terms 
in the Lagrangian are written 

Cx 99 = -\s r {X)F; v F^ - \Sj{X)e^F; v F^ (6) 

where Sr and Si are real and imaginary parts of the gauge kinetic function, respectively. 
In the above, we have taken the S(X) to be universal for all gauge groups, for simplicity. 
Expanding Eq. (6) around the VEV of the X field (X), we find 

£x 99 = -)(s r )f: u f^-U(^\ )5X r F°F^ + ]((^) SsXjFZF^ 



\\ /til \ \ /If 

-\{Si)e^F; v F^ - l(^^y X ^F^F% - \ ^ J SX^F^F^, (7) 

where 5X = X — (X) . It is straightforward to compute the decay width to the gauge boson 
pairs. The result is 

T(X R - 99 ) = V(X, „) = ^Ijfc, (8) 



where d g is a dimensionless constant of order unity defined (in the Planck unit) by 



d 9 = (Kxx*rns„ 



R 



\dX 



(9) 



and Nq is the number of the gauge bosons. Nq = 12 for the minimal supersymmetric 
standard model. In deriving the above result, we have rescaled the field 5X and the gauge 
fields into canonically normalized ones. In Eq. (8), the reduced Planck scale has been 
explicitly written. To give an example of d g , let us consider the Kahler potential of the form 
(2) and the gauge kinetic function S(X) = X. Then d g = 2/y/n, and in fact it is of order 
unity. 

Evaluation of the decay width into the gaugino pairs can be done similarly. We denote 
the gaugino fields by A*") and in two component formalism. The relevant terms of the 
Lagrangian in this case are 



C xss = ^S fi (A')[A<°V"C ( ,A«°l + W^Al-'l - \s,(X)V r [A«V'A W 



+i§^ <o ' m (io) 
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where the covariant derivative is defined as 



P„A« = + \ (Kjd^ - K r d^* j ) A<°> + ■■■. (11) 

Utilizing the equations of motion for the gauginos, one finds that the first line of Eq. (10) 
makes small contributions to the decay amplitude, which are suppressed by the small gaugino 
masses. On the other hand, the contributions from the second line are unsuppressed. In 
fact 




dX \dX I \dX\dX J I \dX* \dX 

, v : v , / dS dF x 

' ^dx^x 6+ \dxdx^ 





*(§*)-(£}•*•»- (§)— . 



where use of Eqs. (4) and (5) is made to obtain the last equality. In the case where mx 3> 
m 3 / 2 , the decay width is simply given 2 



F(X R -> gg) = r(Xj -> gg) = J^d 2 ^. (13) 

128tt 9 M|, V ; 



Notice that it is identical to the decay width to the gauge bosons. When the gravitino mass 
is comparable to the modulus mass, the above result is modified, but remains the same order 
of magnitude. 

The two-body decays of the modulus field into Standard Model fermion pairs as well as 
sfermions can be shown to be suppressed by powers of the masses of final states by using 
their equations of motion [6]. 3 The three-body decays such as quark-quark-gluon will not 
receive this chiral suppression. However they will be suppressed by a s /ir and they will not 
dominate over the decays into the gauge bosons and gauginos. Thus we discard them in 
the subsequent discussion. On the other hand, terms like X^H u Hd in the Kahler potential, 
with H u and H d being the Higgs multiplets in the MSSM, may make sizable contribution to 
the decay width, which can be comparable to the gauge and gaugino final states. Whether 
these terms exist or not are quite model dependent, and we will not consider them here. 

We now examine the modulus decay into the gravitino pair. In the Unitary gauge, the 
relevant interaction terms are in the gravitino bilinear terms in the supergravity Lagrangian 

A/2 = -e^^a u V p ^ - e K ' 2 W*^^ v - e K / 2 W^a^ u , (14) 

where stands for the gravitino in two component formalism, and the covariant derivative 
is given by 



2 The contributions from the F-auxiliary part have not been discussed in Ref. [6]. 

3 For the decays into sfermion pairs, there are also contributions which are suppressed by powers 
of m s /2/mx- These are irrelevant as the gravitino mass is much smaller than the modulus mass. 
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V p i, a = d p ^ + ^c^' - Kj.d^)^ + ■■■. (15) 
Making a field-dependent chiral transformation 

/ w \ ~ 1/4 

Eq. (14) reduces to 

-e^^c^ - ^^(G^ - GrdpP*)^^ - e G ^a^ v + ^a^), (17) 
where G is the total Kahler potential defined by 

G(X,X*) = K(X,X*) + \og\W(X)\ 2 (18) 

In this convention, the gravitino mass is = (e G ^ 2 ). By expanding the Lagrangian (17) 
in terms of 5X and 5X*, one finds the interaction terms 

-\e^{{G x )d p 8X - (G x ,)d p 5X*)^a^ a 

-\m, /2 {{G x )5X + {G x *)5X*){ii^i, v +^,a^i, u ). (19) 

We note that the coupling is governed by (G x ). It is related to the auxiliary field expectation 
value as follows: 

(F x ) = -{{G xx *)- l e G ' 2 G x *) (20) 

Since the VEV of F x is naturally of order m\, 2 /m x , one expects 

(Gx) ~ m 3/2 /m x (21) 

in the Planck unit, which we assume to be the case in the following discussion. We should 
stress that this is indeed the case for the KKLT-type set-up and also for the race-track 
supersymmetry breaking scenario. 

We are now at the position to compute the decay width into the gravitino pair. It is 
evident that the helicity ±1/2 components will give dominant contributions if they are non- 
vanishing, because they will contain enhancement factor of l/m 3 / 2 . To see the point, let us 
take, for simplicity, (G x ) to be real and consider the decay of the real component of SX. 
Then the decay amplitude for a given set of helicity components is written 4 



4 The chiral transformation (16) makes the expression of the matrix element very simple. If we 
use the original Lagrangian (14), then the amplitude contains more than one term. In this case we 
checked that partial cancellation takes place between different terms, arriving at the same result 
we describe here. 
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M(X R - ^3/2^3/2) = i^(Gxx.)- 1/2 (e G/2 Gx)v IM (k')u' i (k), (22) 

where and are gravitino wave functions (in four component formalism). To derive 
this, we have used that the gravitino is a Majorana fermion and the two wave functions are 
related by the Majorana condition = CuJ^ with C being the charge conjugation matrix. 
The wave functions of a massive spin 3/2 field are conveniently expressed as a tensor product 
of a vector and a spinor. For instance, the helicity 1/2 component is written 

u^k; 1/2) = ^(k; 0)u(k; 1/2) + ^(k; l)u(k; -1/2) (23) 

with self-explanatory notation. Here e At (fc, 0) ~ A^/to 3 / 2 at a high-energy limit. The decay 
amplitude into the two helicity 1/2 components of the gravitino is thus expressed as 

M+i/2 = M(X R -> ^3/2(1/2)^3/2(1/2)) 



^(Gxx*)- 1/2 (e G/2 G x )e;(A;'; 0)e*(fc; 0)v(k'; l/2)«(fc; 1/2) 
~ i^Gxx^ieCI'Gx) 1 ^^'; l/2)«(fc; 1/2) 



'V2 



m 3/2 

1 



- d 3/2 m X ( 24 ) 



where <i 3 / 2 is defined as follows: 
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(G^)- 1/2 (e G / 2 G x )=rf 3 / 2 ^. (25) 

m x 

As was seen previously, rf 3 / 2 is a dimensionless constant of order unity. Notice that the l.h.s. 
of the above is the F-auxiliary field of the canonical normalized supermultiplet. The same 
expression is obtained for the decay into the helicity —1/2 components. On the other hand, 
the decays into the helicity ±3/2 components are suppressed by powers of m 3 / 2 /m X - Thus 
the decay width is computed to be 



m 
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UX B -+ = (26) 

at the limit m 3 / 2 <C mx- Here the reduced Planck scale has been written explicitly in the 
final expression of the decay width. A computation can also be performed for the imaginary 
part of the 8X, with the same result as Eq. (26). 

Thus the modulus field dominantly decays into the gauge bosons and the gauginos with 
the total decay width 



tot 



F(X - all) ~ F(X - gg) + Y{X - gg) 

Ng 3 ( Nq\ 2 m\ 

a g i\j2 \ zi ) 



64tt 9 16?r V 12 J 9 

and the branching ratio of the decay to the gravitino pair 



7 



1 / Nr\ _1 dn/n 

B 3/2 = Br(X - ^3/2) = Va ) (28) 

With d g and rf 3 / 2 being the constants of order unity, we find that the branching ratio to the 
gravitinos is of order 1CT 2 . 

The production of the gravitinos at the modulus decay has striking impacts on the 
cosmology. Here we consider the situation where the coherent oscillation of the modulus 
field will dominate the energy density of the universe after primordial inflation, releasing 
huge entropy and reheating the universe at the decay. The reheating temperature at the 
modulus decay is estimated by equating the total decay width to the expansion rate of the 
universe at the reheating: 

rp ( 90 

J-R 




where g*{Tn) is the effective degrees of freedom of the radiation at the reheating. The 
gravitino yield produced by the modulus decay, which is defined by the ratio of the gravitino 
number density n 3 / 2 relative to the entropy density s, can easily be evaluated as 

3„ T R 



-B 



-3/2 

2 m x 



= 0.-3 x io- Bs/2 feewv" 4 (^) U2 dl (-^-f 

61 I 10 J V 12 7 9 V10 5 GeV/ 



(30) 



A constraint on the gravitino yield comes from the big-bang nucleosynthesis (BBN). As 
is well-known, the success of the BBN would be threatened by the electromagnetic showers 
(as well as the hadronic showers) produced at the gravitino decay. It is termed the gravitino 
problem. The gravitinos are produced via scattering processes in the thermal bath after 
primordial inflation epoch. In the situation we are considering, the gravitinos in this origin 
are diluted by the entropy production at the modulus decay. However they are regenerated 
directly by the modulus decay. The requirement that the gravitino decay products should 
not spoil the BBN severely constrains the gravitino abundance. Recent analyses [23-25] 
show that, for the gravitino mass in the range m 3 /2 — 10 3 — 10 4 GeV, the constraints from 
D/H and 6 Li are the severest, leading to Y3/2 ^ 10 -16 when the hadronic branching ratio of 
the gravitino decay B h is 1, and Y 3 / 2 ~ 10~ 13 when B h = 10~ 3 . Comparing these numbers 
with the gravitino yield obtained at the modulus decay (30), one sees that the latter exceeds 
by several orders of magnitude. For lighter gravitino m 3 / 2 ~ 10 2 — 10 3 GeV, 3 He/D also 
plays a role. The constraint in this range remains very severe, roughly speaking at the level 
^3/2 ~ 10 -16 . It becomes somewhat weaker for heavier gravitino, especially for the case 
Bh = 10 -3 . However the yield (30) still exceeds the constraint from 4 He abundance. Finally 
the constraint disappears when the mass is above 10 5 GeV, corresponding to the life-time 
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shorter about 1CT 2 sec. Thus we conclude that the (unstable) gravitino whose mass is less 
than 10 5 GeV is excluded by the BBN constraint. 

Though we already obtained a very severe bound on the gravitino mass, this is not 
the end of the story. When the gravitino is unstable, it decays to lighter superparticles, 
i.e. -R-parity odd particles. Under the assumption of imparity conservation, the lightest 
superparticle (LSP) is absolutely stable. A plausible candidate for the LSP in the MSSM is 
the lightest in the neutralino sector, that is a linear combination of the neutral gauginos and 
higgsinos. The neutralino LSP, if stable, will contribute as (a part of) the cold dark matter 
whose abundance is bounded from cosmological observations. The second question we would 
like to address is therefore whether the neutralino abundance produced by the gravitino 
decay will not exceed the upperbound of the dark matter inferred by the observations. We 
should note here that the problem of the over-abundance of the neutralinos produced by the 
moduli decay was addressed in Refs. [3,5]. 

The gravitinos produced at the modulus decay do not interact with others. Thus the 
yield of the gravitinos does not change until they decay. The gravitino decay width into the 
MSSM particles is (see for instance Ref. [26]) 



13/2 " 384^' ( } 



when all MSSM (super)particles are included in the final states and their masses are ne- 
glected, which is justified for the gravitino much heavier than the MSSM (super)particles. 
It follows from Eq. (31) that the temperature of the universe at the gravitino decay is 



90 



.1 



1/4 



Numerically it reads 



r3/2 , 7 , xl0 -3(^M)- I/4 (_^_)- Ge V ,33, 

The gravitino decays into lighter superparticles, followed by cascade decays to the neu- 
tralino LSP. The neutralino LSPs produced this way are so abundant that they annihilate 
with each other. The annihilation process terminates when the annihilation rate reduces to 
the expansion rate of the universe at the gravitino decay if (T 3 / 2 ) [3] 

(Oann^rcl)™* - H ( T 3/2), (34) 

where <r ann is the annihilation cross section of the two neutralino LSPs, t> rc i their relative 
velocity, (• • •) represents the average over the LSP momentum distribution, and n x is the 
number density of the neutralino LSPs. The above argument derives an estimate for the 
neutralino abundance: 



H(T 3/2 ) 



T 3/2 (Cann^rel)s 



90 \ 



1/2 



T 3/2 



n 2 g*(T 3/2 ) ) (o- ann w re i)T3 /2 M P1 



(35) 
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The yield remains constant until today. A more sophisticated evaluation requires to solve 
Boltzmann equations numerically or analytically [27], but the estimate given above is suffi- 
cient for the purpose of the present paper. 

The annihilation cross section depends on the LSP component as well as the superparticle 
mass spectrum. To maximize the annihilation effects in the MSSM, let us consider the wino 
LSP, the neutral component of the SU(2)l gauginos. Assuming that it is heavier than the 
VT-boson, the dominant mode of the wino annihilation is into W-boson pair via charged 
wino exchange. The annihilation cross section is computed to be [6] 

2nm^ (2 — xw> 

where g 2 is the SU{2) L gauge coupling constant, m x the wino mass, and x w = rri^/m 2 ,. 
Here possible co-annihilation effect has not been taken into account [28]. 

With (35) and (36), it is straightforward to compute the wino relic abundance today. 
The result is conveniently expressed by the relic mass density relative to the entropy density: 

~ o 43 x 10- 9 GcV ( 2 -^ 2 ( 9*( T ^ \ ( m x \ 3 ( m ^ \ ^ m 
s " U - 4dxiU GeV (i _ Xw) s/2 [ 10 ) UooGevJ U 5 GeW ' {67) 

or in terms of the density parameter fl x which is defined by the ratio of the LSP mass density 
to the critical mass density of the universe 



xtl ~ U (1 - xwfl* \ 10 ) UoOGevMlOSGevJ 



(38) 



with h ~ 0.72 being the Hubble constant in units of 100 km/s/Mpc. 

In Fig. 1, the constant contours of the density parameter fl x h 2 are drawn in the m x -m 3 / 2 
plane. The real lines represent Q x h 2 =0.01, 0.1, and 1. Given a LSP mass, the density 
parameter decreases as the gravitino mass increases. In the same figure, we also show the 
contour of fi x h 2 = 0.13, roughly corresponding to the 95%CL upperbound of the cold dark 
matter abundance from the cosmological observations [29] . Thus in order to avoid too much 
abundance of the neutralino LSPs, the region below this line should be excluded. We find 
that for the wino mass in the weak scale the gravitino should be heavier than (9(10 6 ) GeV. 
On the other hand, around the line of Q x h 2 = 0.1, the wino-like LSPs produced by the 
gravitino can constitute the dark matter of the universe. Here we should caution that the 
estimate of the relic abundance given in this paper is rather rough, which may contain an 
error of factor 2 or so. 

As we mentioned earlier, the wino case gives the largest annihilation cross section and 
thus the weakest constraint on the gravitino mass. A similar but slightly severer bound will 
be obtained for the higgsino LSP case. On the other hand, the bino LSPs will remain too 
much with such a low T 3 / 2 [3,5]. We note that when the gravitino mass becomes heavier than 
about 10 7 GeV, the produced LSPs will get thermalized and the conventional computation 
of the relic abundance can apply. 

Here we would like to briefly discuss what would happen if the decay into the gravitinos 
is negligibly small. This is the case when the VEV of auxiliary field Fx of X is (accidentally) 
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small. In this case the superparticles are produced at the modulus decay with the branching 
ratio of 0.5. A similar argument given for the gravitino decay can apply except that the 
gravitino decay temperature should be replaced with the reheating temperature of the mod- 
ulus field. Then we obtain quite a similar lower bound on the modulus mass to avoid the 
over- abundance of the neutralino LSPs [3,5]. However there are some substantial differences 
in the two cases. The modulus mass of the order 10 6 or 10 7 GeV may be acceptable from 
model building point of view when one considers a KKLT-type model or a racetrack model. 
There are models in which the modulus mass is separated from the supersymmetry breaking 
scale. A simple model to realize this situation was given in Ref. [30]. On the other hand, 
the gravitino mass is directly related to the supersymmetry breaking scale. Thus one has to 
elaborate model building to accord with a very heavy gravitino. We will come back to this 
point shortly. 

Let us summarize what we have obtained in this paper. We have considered the moduli 
decays into various two body final states and discussed the impacts on the cosmology. It 
turns out that the total decay width is of order m^-/167rMpj as far as the coupling of the 
modulus field to the gauge multiplets is not suppressed in the gauge kinetic function. The 
main decay modes in this case are the decays into the gauge boson pairs as well as those into 
the gauginos. The relevant coupling to the latter decay modes emerge when the auxiliary 
component of the modulus field is integrated. The most important is the decay into the 
gravitino pair. We have shown that the coupling of the modulus field to the helicity ±1/2 
components of the gravitino is proportional to the VEV of the modulus auxiliary field, i.e. 
the supersymmetry breaking of the modulus field. It is known that this VEV is of the order 
m\i 2 /mx in the Planck unit for the field with a Planck scale expectation value. Then the 
decay width into the gravitino pair follows a naive dimensional counting and its branching 
ratio is at (a few) % level. This large branching ratio has striking impacts on the cosmology 
associated with the gravitino. Firstly the daughter particles produced by the gravitino 
decay would upset the BBN. What is worse, the neutralino LSPs produced by the gravitino 
decays would exceed the abundances of the cold dark matter inferred by the cosmological 
observations. These constraints require the mass of the gravitino heavier than at least about 
10 6 GeV for the wino-like LSP case. For other types of neutralino LSPs, the gravitino mass 
bound may reach 10 7 GeV. 

Such a heavy gravitino will cause serious difficulty to model building of low-energy super- 
symmetry. In particular, superconformal anomaly mediation [31] will make contributions 
to soft masses. They are typically proportional to the VEV of the auxiliary component 
of the superconformal compensator F^, suppressed by one-loop factor of 1/IQtc 2 . Naively 
the value of F^ is comparable to the gravitino mass. If this is the case, with the gravitino 
mass of 10 6 GeV or higher, the resulting soft masses would be far above the electroweak 
scale, diminishing the very motivation of low-energy supersymmetry. One may speculate 
that compactification with warped metric may help, where the MSSM sector is localized on 
a brane with a warp factor [32] (see also Ref. [33]). A modest warp factor of 10 _1 — 10~ 2 
will be enough. Then the VEV of F^ at the location of the MSSM brane will receive the 
suppression from the warp factor compared to the gravitino mass. 
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FIGURES 
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FIG. 1. Constant contours of the density parameter Q x h 2 are drawn in the m x -m^/2 plane. 
Three real lines represent ^ x /i 2 =0.01, 0.1, and 1, from the above. Also shown by a dashed line is 
the contour of Q x h 2 = 0.13, approximately corresponding to the 95%CL upperbound of the cold 
dark matter abundance from the cosmological observations. 
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